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ABSTRACT 

We report rest-frame optical Keck/NIRSPEC spectroscopy of the bright lensed galaxy RCSGA 
032727-132609 at z=1.7037. From precise measurements of the nebular lines, we infer a number of 
physical properties: redshift, extinction, star formation rate, ionization parameter, electron density, 
electron temperature, oxygen abundance, and N/0, Ne/0, and Ar/0 abundance ratios. The limit on 
[O HI] 4363 A tightly constrains the oxygen abundance via the "direct" or Te method, for the first 
time in an average-metallicity galaxy at Z'^2. We compare this result to several standard "bright-line" 
O abundance diagnostics, thereby testing these empirically-calibrated diagnostics in situ. Finally, we 
explore the positions of lensed and unlensed galaxies in standard diagnostic diagrams, and explore 
the diversity of ionization conditions and mass-metallicity ratios at z=2. 
Subject headings: galaxies: high-redshift — galaxies: evolution — gravitational lensing 
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1. INTRODUCTION 

Our knowledge of the Universe's star formation his- 
tory has advanced remarkab ly in the past doze n years . 
First rest-UV photometry (jMadau et all 119961 I1998D . 
then 24 ^m Spitzer photometry for < ^ < 2 galaxies 
(|Caputi et al.ll2007l: ILe Floc'h et "all 120051) and rest-UV 
photometry at high er redshift (c.f. lOuchi 
iBouwens et al.l [20101) . has shown us that the SFR rises 
steeply from z=0 to z=l, has a broad plateau from z=l to 
z=4-5, and falls off (with debated slope) out to reioniza- 
tion. Thus, the question, "What is the star formation 
history of the Universe?" has been answered reason- 
ably, and the focus has shifted toward, "Which galax- 
ies formed their stars when, and why, and how did that 
star-formation change those galaxies and their environ- 
ments?" 

While we have identified the galaxies that formed 
most of the Universe's stars, we have much to learn 
about how that process occurred. One key question is, 
what were the physical conditions inside these galaxies- 
metallicity, abundance, extinction, stellar effective tem- 
perature, electron temperature and density — compared 
to star-forming galaxies today? How did these physical 
conditions evolve through episodes of star formation and 
gas accretion? 

Until the era of extremely large telescopes arrives, our 
best chance to address this question is to study galax- 
ies that are highly magnified by gravitational lensing. In 
such rare cases, magnification factors of 20-30 make di- 
agnostic spectroscopy possible with current telescopes. 
Such work was pioneered in the g alaxy MS151 2 -cB58 
(|Yee et al.lll996l : hereafter cB58) bv lPettini et"all (pOOOl ) 



and iTeplitz et al.l (|2000D . and can now be extended to 
a larger sample thanks to disco veries of bright l ensed 
galaxies (e.g., lA Ham et al. 2007 : iBelokurov et al] 120071 : 
iSmail et al.ll2007i:iRigbv et al.,i2008HKoester et aLlboiOl ) 
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Among these new discoveries is RCS GA 032727- 
13260 9, hereafter RCS0327, at z=1.7 (jWuvts et al.l 
l2010f ). discovered from the second Red Sequence Cluster 
Survey (RCS2; Gilbank et al. in prep.). With an inte- 
grated g-band magnitude of 19.15, we believe RCS0327 
to be the brightest high-redshift lensed galaxy yet found. 
iWuvts et al.r()2010[ ) present the discovery, spectroscopic 
confirmation, deep-follow-up imaging, spectral energy 
distribution modeling, and preliminary lensing analysis. 
In this paper, using 1.3 hr of integration with NIRSPEC 
on Keck, we determine with unprecedented precision the 
physical conditions of star formation in this hopefully- 
typical galaxy at the crucial epoch of z~2. 

2. METHODS 

RCS0327 was observed on 04 Feb. 2010 UT with the 
NIRSPEC spectrograph (McLean et al. 1998) on the 
Keck II telescope. The weather was clear, and the seeing 
was measured as 0.85" and 0.45" when the telescope was 
focused during the night. We used the low-resolution 
mode and the 0.76" x 42" longslit. We targeted the 
brightest ^ 10" of the arc, at a position angle of 134°, as 
shown in Figure [TJ The target was acquired by offsetting 
from nearby stars on the near-IR slit-viewing camera, 
target acquisition was verified by direct imaging on this 
camera. The target was nodded along the slit in an AB 
pattern, with exposures of 600 s per nod. Table 1 sum- 
marizes the filters and exposure times. The AOV star HD 
23683 was observed every hour as a telluric standard. 

We reduced the spectra using the nirspec-reduce pack- 
age written by one of us (G. D. Becker), which uses lamp 
exposures to fiatten the data, the sky lines to wavelength 
calibra te, and optimally fits and subtracts the sky fol- 
lowing |KeIio3 (|2003). For each frame, we measured the 
spatial profile of the lensed arc by fitting the brightest 
emission lines, then used this spatial profile to optimally 
extract the spectrum. 

The arc is extended over 38", roughly 10" of which 
was captured by the NIRSPEC slit. In a subsequent 
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filter 

NIRSPEC-1 
NIRSPEC-3 
NIRSPEC-6 



Table 1 

Observation Log 

t(s) wavelength range(/im) 



Table 2 

Measured line fluxes for RCS0327 



1200 
2400 
1200 



0.948-1.16 

1.14-1.43 

1.76-2.19 



paper, we will analyze the spatial variation of physical 
conditions across the arc. In this paper, we consider the 
integrated spectrum. 

Each extracted spectrum was corrected for telluric 
absorption and flu xed using the tool xtellcor_general 
(|Vacca et al.l [200f). using the closest-in-time observa- 
tion of the AOV standard star. The flux level is thus 
appropriate for the fraction of the galaxy inside the slit, 
not for the whole galaxy. In we estimate the fac- 
tor by which our NIRSPEC fluxes should be scaled to 
represent the whole galaxy. 

Using a telluric standard to flux the spectra provides 
excellent relative fluxing within a given fllter, but be- 
cause observations were taken in three separate fllters, 
there can be offsets between fllters due to differential slit 
losses, due for example to changes in seeing or pointing. 
In addition, the lines in fllter NIRSPEC-6 are observed 
right at the edge of an atmospheric transmission window, 
and thus may suffer especially high telluric variability. 
We address these issues in i i3.2l 




Figure 1. Finder chart for RCS0327. This J-band image from 
IWuvts et an 1 120101 illustrates how we positioned the NIRSPEC 
slit, which is 42 X 0.76". N is up and E is left. For simplicity 
we show the slit centered on the source; in fact it the source was 
placed on the left half and then the right half of the slit ( "an AB 
pattern"), with the nods separated by 10—15". 

For each filter, we combined the individual flux- 
calibrated ID spectra with a weighted average, producing 
one fluxed spectrum for each of the three filters. 

We fit line fiuxes as follows. We fit each isolated 
line with a Gaussian to measure the line flux, using 
Levenberg-Marquardt fltting. The continuum used was 
the mean flux in an adjacent spectral region, chosen by 
hand for each line. For partially-overlapping lines, we 
simultaneously fit multiple Gaussians using A. Marble's 
implementation of the IDL Levenberg-Marquardt least- 
squares fitting code MPFITFUN (,Markwardt.2009.) . We 



line 



flux dflux 



Filter Nl 

[O II] 3727 
[Ne III] 3869 
H8-I-Hel 3890 
He 
H<5 

Filter N3 

H7 

[O III] 4363 
[Ar IV] 4741 
H/3 

[O III] 4959 
[O III] 5007 
Filter N6 

[N II] 6548 
Ha 

[N II] 6583 
[S II] 6716 
[Ar III] 7136 



1.00808 

1.046273 

1.05195 

1.0729 

1.1093 

1.17385 

1.1800 

1.2825 

1.31473 

1.34111 

1.35409 

1.77115 
1.77488 
1.78002 
1.81663 
1.9296 



72.0 

7.8 

6.6 

<3.8 

7.1 

13.5 
<1.5 
<3.0 
32.4 
49.3 
159 

3.85 

116 

7.4 

6.7 

2.8 



0.2 

1.3 

0.4 

limit 

1.7 

2 

limit 

limit 

1.1 

1.6 

1.4 

0.8 
1.2 
1.2 
0.3 
0.3 



Note. — Measured line fluxes. 
Columns are line ID, observed 
wavelength, lineflux in units of 
10~^^ erg s~^ cm"'^, and uncertainty in 
the lineflux. 

report line fluxes in Tabled) 
We assume a cosmology of V, 

-1 C^1„ 



0.3, riA = 0.7, and 
Hq = 70 km s ~^ Mpc~^. S olar abundances are taken 
from Table 1 of lAsplund eFal . (2009). The initial mass 
function (IMF) is iChabrieii (|2003i) unless noted. 



3. RESULTS 

Figured] plots the spectra. Several of the nebular emis- 
sion lines are detected at very high signal-to-noise ratio 
(SNR); for example, we detect ~50,000 net counts in 
the Hz? line. Continuum is detected in filters NIRSPEC- 
1 and NIRSPEC-3 at SNR < 1 per pixel. Continuum 
is not confidently detected near the lines of interest in 
NIRSPEC-6, which is unsurprising given the low atmo- 
spheric transmission at those wavelengths. 

3.1. Extinction 

Based on SED fitting, IWuvts et all (|2010D report low 
extinction values of E(B-V) = 0.03 to 0.11, depending on 
metallicity and spatial position. Given the potential de- 
generacies when fitting SEDs between extinction and age 
and star formation rate, a spectroscopic measurement of 
the extinction is desired. 

The a, /3, 7, and S transitions of the Balmer series are 
detected; He is formally undetected. The c orrection for 
stellar absorption, ^2A (jMcCall et al.ll985l) , is negligible 
at these extremely high equivalent widths: E^ = 430 ± 
70 A and 100 ± 20 A for and H7. 

Ha appears in fiflcr NIRSPEC-6; H7 and H/3 in 
NIRSPEC-3; and HS and He in NIRSPEC-1. A line is 
detected at the position of H8, but is too strong to be 
H8 alone. We suspect it is a blend of H8 with tw o He I 
lines, as in seen in the Orion Nebula tOstcrbroc k et al.l 
[1991) . 

Because they appear in the same filter, the H/3/H7 and 
H(5/He line ratios are free from relative fluxing errors. 
We measure H/3/H7 = 2.39 ± 0.37, which for Case B 
recombination at T= 10^ K and Ue = 100 cm^'^ yields 
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E(B-V) = 0.23 ± 0.2^1 The US/He ratio of > 1.87 does 
not constrain the extinction given measurement errors. 
The best spectroscopic measure of extinction therefore 
comes from the H/3/H7 ratio, which for R=3.1, yields 
A„ = 0.7 ±0.7. 

The high end of this range is inconsistent with the 
blue color of the ar c, as quantified by the SED fitting of 
iWuvts et al] ()2010[) . Therefore, in the rest of the paper 
we will report results assuming first Ay — 0.0 and then 
Ay = 0.7. This considerable uncertainty in the extinction 
will limit how well we can measure spectral diagnostics 
that span a large range of wavelength, for example R23. 
A more precise measurement of the extinction awaits si- 
multaneous observations of Ha and H/3, or a longer inte- 
gration on H7 and H/3. 

3.2. Tweaking the flux calibration 

As discussed in ij2l there may be fluxing offsets among 
the three filters. We address this issue using the Balmer 
series. Assuming an extinction of A„ = 0.7 as measured 
from the H/3/H7 ratio in NIRSPEC-3, we tweak the rela- 
tive scalings of the NIRSPEC-1 and NIRSPEC-6 spectra 
to bring Ha/H/3 and H(5/H/3 to the appropriate ratios for 
this extinction. This increases the flux in NIRSPEC-1 
by a factor of 1.15 relative to NIRSPEC-3, and decreases 
the flux in NIRSPEC-6 by a factor of 0.61 relative to 
NIRSPEC-3. These are reasonable corrections given the 
slit loss variation expected from seeing changes and at- 
mospheric transparency issues in NIRSPEC-6. Fluxes 
reported in this paper reflect this scaling. 

In TableHwe report the line fluxes in filters NIRSPEC- 
1,-3, and -6. Line flux ratios within a filter should be 
precise, limited by the uncertainty in line-fitting. When 
calculating a line flux ratio for lines spanning different 
filters, one should include a 10% relative fiuxing uncer- 
tainty. 

3.3. Fraction of the galaxy covered 

To estimate the fraction of the galaxy covered by the 
slit, we use the PA NIC /Magellan images in J and Kg of 
IWuvts et all (l20Tol) . which had seeing of FWHM=0.57" 
and 0.51". We estimate that 37 ± 5% of the total arc 
light falls into a slit positioned as in Figure [H neglecting 
slit losses. We compute slit losses for seeing of 0.8-0.9" 
and pointing errors of 0-0.25". We conclude that 32% ± 
4% (random) ±1% (systematic) of the continuum light 
should have been captured by our NIRSPEC slit. 

We now estimate this quantity in an alternate way. We 
measure the continuum adjacent to the 3737 A and H/3 
lines as 6.1 x 10~^* and 6.0 x 10~^^ erg s~^ cm~^ 
(using a robust average), with measurement errors of 10- 
20%. These continuua levels are 21% of the continuum 
in tensity in the spect ral energy distribution in Figure 7 
of IWuvts et"all (|2010[ ). which is the best fit to the u, B, 
g, r, I, z, J, H, Ks photometry of the whole arc. 

These two methods disagree by 38%. We suspect that 
the continuum method is as fault, as the continuum is not 
well detected, and small DC offsets might plausibly be 
introduced in sky subtraction. (Indeed, this uncertainty 
in the continuum level is why we chose to fiux using AOV 
stars rather than multiplying line equivalent widths by 



f^ from broadband photometry as in iTeplitz et ai][2000l ) 
Therefore, we adopt the first method, and will multiply 
our NIRSPEC fiuxes by the reciprocal, 3.12 ± 0.4, to 
convert to the expected line flux over the whole arc. 

3.4. Average Magnification 

IWuvts et all ([20lol) measure an average magnification 
of It = 17.2 ± 1.4 for RCS0327. Our cu rrent lensing 
mode l, based on ground-based photometry (jWuvts et al.l 
'2010) , has insufficient spatial resolution to determine the 
magnification of each small knot. Thus, at present we 
cannot quote a reliable magnification for the portion of 
the arc that falls in our NIRSPEC slit. It is likely that 
some of the bright knots are much more highly magni- 
fied than the arc on average. For the time being we 
will adopt the average magnification, with these caveats. 
Fortunately, except for the star formation rate derived 
from the Balmer lines, none of the results in this paper 
depend on the magnification within the NIRSPEC slit. 

3.5. Redshift 

We fit the mean emission line redshift as 1.7037±0.0001 
using the NIST linelistH, excluding [O II] 3727 be- 
cause it is a partially resolved doublet. This compares 
perfectly with the redshift of rest-UV emission lines 
in our (not yet published) MagE/Magellan spectrum: 
1.70369 ± 0.00006. The MagE spectrum also yields a 
redshift for absorption lines in the interstellar medium: 
Sill 1260, Si II 1808, and Al II 1670, recommended to us 
by C. Tremonti as isolated ISM lines, yield an ISM red- 
shift of 1.702671 ±0.0003. Thus, the MagE spectra show 
that the interstellar medium is blueshifted from the neb- 
ular lines by 110 ± 30 km s~^. Our measured redshifts 
are 3 and 2a higher than the 1.7009 ± 0.0008 reported 
by Wuyts et al. (2010) using a combination of the [O II] 
3727 and C HI] 1909 emission lines and the Fe II and 
Mg II absorption lines. 

3.6. Velocity width 

We extract arc lamp spectra using the same aper- 
tures as the science frames, and combine with the same 
weighted average technique as for the science frames. For 
each filter, we fit the instrumental linewidth versus wave- 
length relation with a linear fit, then interpolate the in- 
strumental linewidth at the wavelength of each astro- 
nomical line of interest. For example, at the observed 
wavelength of Ha, the arc fines have FWHM=13.6A 
compared to FWHM=15.3A observed for Ha; from this 
we infer cr = 50 ± 2 km s'^ For H7, H/3, A4959, and 
A5007, we measure a = 76, 59, 49, 23 km s~\ With Ha, 
this gives an average nebular line velocity dispersion of 
(7 = 51 km s~^ with an error in the mean of 9 km s~^. 

We fit the [O II] 3726, 3729 doublet with two Gaus- 
sians, fixing the line centers using the NIST vacuum 
wavelengths and redshift from i i3.5l and forcing the two 
linewidths to be the same, while varying the continuum, 
line amplitudes, and the linewidth. This fitting finds that 
cr = 51 ± 1 km s-i for [O II]. 

Because the lensing morphology of RCS0327 is com- 
plex, we cannot yet compute a meaningful effective ra- 
dius, as needed to determine a dynamical mass from the 



Using the extinction law of lCardelli et all l|1989l) . 



''' http://www.pa.uky.edu/ peter/atomic/ 



4 



Rigby et al. 



velocity width. Upcoming HST observations should clar- 
ify this matter. 

3.7. Star formation rate 

iWuvts et ail (|2010( ) estimated the SFR by fitting the 
broad-band optical and near-IR photometry, finding a 
SFR constraint of < 77Mq yr~^ for 40% solar metahic- 
ity. (For solar metallicity the SFR is lower). 

We also estimate the star formation rate from mid- 
IR photometry. We measure the main arc and counter- 
image separately, finding MlPS/Spitzer 24 /um flux den- 
sities of 1040 ± 153 /zJy and 43 ± 11 /zJy respectively. 
Scaling by the average magnifications for the main 
arc and counter-im age from the current lensing model 
(jWuvts et al.ll2010D yields de-lensed 24 /im flux densi- 
ties of 60 ± 10 /LtJy for the main arc versus 21 ± 6 /iJy 
for the counter-image. The de-lensed 24 yitm flux density 
of the main arc corresponds to a star form a.tion rate of 
106 lb 3OM0 yr^^ using the prescription of iRieke et al.l 
()2009() . The equivalent value from the counter-arc is 
18±8 M© yr-i. 

The significant difference between the de-lensed flux 
densities suggests that the source is strongly non- 
uniform; resolving this requires a more precise lensing 
model and a robust comparison of the observed image- 
plane data (both spectra and photometry) in the source 
plane. Until HST imaging of this source is available, the 
best comparison to the spectra presented here is the main 
arc, as the spectra sample a portion of that image. 

Thirdly, we estimate the star formation rate from the 
NIRSPEC spectra, using the H/3 line flux within the 
aperture, assuming the extinction f rom H/3/H7, and the 
S FR(Ha) convers ion of iKennicuttI ([l998. i. modified for 
alChabrieii (|2003l) IMF following Eqn. 10 of IRieke et all 
(I2OOI I 

L{Ha) = f{Hl3){Ha/Hp)ATrdl/u (1) 
= (3.2 ± 0.1) X 10^1^(3.6 ± 0.8)47rdi/M (2) 
= (2.25 ± 0.5) X W'^/u (3) 
= (5.9± 1.3) X lO^Vw^O (4) 

SFR{Ha) = 0.66 x 7.9 x lO^^'^erg (5) 
= (1170 ±27OM0 ?;r"i)/it (6) 

A rough method of computing the total star forma- 
tion rate is to divid e by the average magnification from 
IWuvts et al.l ()2010D and multiply by the fraction of the 
light in the image plane captured by the NIRSPEC slit 
(from g331). This yields SFR = 21O±6OM0 yr^i, which 
is extremely high. 

The star formation rate derived from the Balmer lines 
is considerably higher than inferred from the broadband 
photometry or the 24 /im flux. The cause is not measure- 
ment error nor fluxing error. (We have cross-checked the 
fluxing in multiple ways.) The Balmer lines are simply 
incredibly bright in these knots of RCS0327, with H/? for 
example eight times brighter than in cB58. 

By spectroscopically targeting the brightest knots in 
RCS0327, we may well be biased toward the regions with 
highest surface brightness and/or highest magnification. 



The current lensing map (jWuvts et all I2O1O0 lacks the 
high spatial resolution required to determine the magni- 
fication of each knot. However, it does suggest that the 
magnification of this part of the arc is significantly higher 
than other parts of the arc. Thus, the rough method of 
computing the total SFR is probably wrong because it 
uses the average magnification of the whole arc, rather 
than the magnification of each pixel within the NIRSPEC 
slit. 

A simpler method would be to measure the emission 
line flux over the entire arc, or at least all of the counter- 
arc, which is more compact and thus easier to map. 
Pending narrow-band HST imaging of H/3 should pro- 
vide this coverage for RCS0327. 

The three different methods give quite different an- 
swers for the star formation rate. Only the brightest 
portion of the main arc was observed with Nirspec; this 
portion of the arc is possibly the brightest because of a 
conflation of a bright region within the source with a re- 
gion of larger-than-average magnification, from which we 
would expect the spectral star formation rate to be larger 
than that derived from the average main arc 24 //m fiux 
density, as observed. A robust treatment of this compar- 
ison awaits a refined lensing model and reconstruction of 
source in the source plane. 

Imaging is also biased toward regions of high surface 
brightness, and over-represents certain portions of the 
galaxy, but it does capture the entire image plane, and 
thus should deliver accurate quantities like average mag- 
nification, total magnitudes, and colors. 

3.8. Ionization parameter 

Figure 1 of iKewlev fc Dopital (|2002l ) iUustrates the 
A5007/A3727 fiux ratio as a diagnostic of the ionization 
parameter. For an O abundance of 20-40% solar (on the 
Asplund system), eqn. 12 of IKewlev fc Dopital (|2002j3 
yields an ionization parameter of log U = —2.73 to —2.85 
for Ay = 0.7. For = 0, log U is higher by 0.1 dex. The 
uncertainty due to extinction dominates over the flux un- 
certainty. 

3.9. Electron Density 

The two-component fit to the [O II] doublet in §3.6! 
found a line flux ratio of f(3726/3729) = 0.893 ± 0.024, 
where this uncertainty is dominated by the uncertainty in 
deblending the dou blet. This ratio is density-de pendent 
(c.f. Figure 5.8 of lOsterbrock fc Ferlandll200l ) Using 
stsdas. analysis. temden in IRAFI3 our measurement cor- 
responds to a tight measurement of the density: — 
235tll at Te = 10^ K, and = 252^1° cm^^ at 
1.2 x lO'* K. 

The ratio of the C HI 1907/1909 d oublet also con- 
strains the density (jRubin et al.l |2004 their Figure 2.) 
Fitting the doublet in our Mage/Magellan spectra (Rigby 
et al. in prep.) in the same way as the 3727 doublet, we 
measure a flux ratio of f(1907/1909) = 2.16 ± 0.3. This 
ratio is unphysical by 1.7(7, as the zero-density limit is 

* Which uses the Anders &: Grevessel l| 19891) abundances; we con- 
vert to Asplund et al. (2009). 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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1.65 for pure C13, and 1.51 for pure C12 ([Rubin et al.l 
I2004D . While deeper data should provide a better mea- 
surement of this line ratio, the current measurement of 
the C III] ratio is consistent with the low density regime 
of this diagnostic of Ue < 10^'^ cm~^. 

The doublet ratio of [S II] 6716, 6731 A also 
constrains the electron density (c.f. Fig. 5.8 of 
lOsterbrock fc Ferlandl 120061 ) Unfortunately, the redder 
line is mostly lost to a skyline. Fitting in the same way 
as for the 3727 doublet, we measure a doublet ratio of 
2.4 ± 0.4, which is unphysical by 2.4tT (the zero density 
limit is 1.43). This unphysical result is likely due to the 
skyline contamination; we conclude that in this case we 
cannot use the [S II] ratio to constrain the electron den- 
sity. 

To summarize, the [O II] 3727 doublet ratio provides a 
tight density constraint, which is consistent with the low 
density regime indicated by the C III] 1909 doublet. The 
[S II] doublet is contaminated by a skyline and provides 
no density constraint. 

3.10. Electron Temperature 

The ratio [O III] (A5007-FA4959)/A4363 constrains the 
electron teniperatu re with almost no rig dependence (c.f 
Izotov et al.l l200l Fi gure 5.1 of lOsterbrock fc Ferlandl 
2006D. The ratio in RCS0327 is > 139 for A „ = 0.0 , and 
> 121 for A„ = 0.7. Following llSotov et al.l 1)20061 ) . this 
corresponds to < 1.14 x 10'' K, < 1.20 x 10"* K, and 
< 1.26 X lO'' K for A„ = 0.0, 0.7, and 1.4. 

3.11. Oxygen abundance 

We constrain the oxyge n abundance via t he "direct" 
or "Te" method following llzotov et al.l ()2006t ). using the 
non-detection of [O III] 4363 A to constrain Te, [O 
II] 3727 A and H/3 to constrain {0+/H+) and [O III] 
4959, 5007, and H/3 to constrain {0^+/H^+). Since He 
II 4686 is not detected, we can ignore the contribution 
of 0^+. The result is 12 + \og{0/H) > 8.21 for Ay = 0.0 
and > 8.14 for Ay = 0.7; these are 0.4 8 and 0.55 dex 
lower than the solar value of 8.69 ± 0.05 ([Asplund et al.l 
I2OOI . 

We compare this "direct" lower limit on the oxy- 
gen abundance to results from the bright line diagnos- 
tics. We follow the calibrations for each diagnostic, 
then remove the r elativ e offsets via the conversions of 
IKewlev fc EllisonI (|2008D . so th at all diagnostics are o n 
the system of the N2 index of IPettini fc Pagell (|2004[ ). 
The bright line results are as follows: 

• The N2 index, \ogi[N 1 1] / H a) = -1.19 ± 0.07, 
yields 12 + \og{0/H) = 8.20 ± 0. 04 by the third- 
order polynomial calibration of IPettini fc Pagell 
()2004D . and 8.22 ± 0.04 by their linear fit. Each 
of these two calibrations was rep orted to have a 
la sp read against Te of 0.18 dex ([Pettini fc Pagell 
120041 ). Reddening is irrelevant. 

• The N2 index calibration of iDenicolo et al.l (|2002D 
yields 12 + log(Ci/iJ) = 8.20 ± 0.05. 

• The 03N2 index, \og[{^007 / H fi) / {N 1 1 / H a)] = 
1.88 ± 0.07, yields 12 + \og(0/H) = 8.16 ± 0.02 
by the calibration of IPettini fc Pagell (|2004l) . which 



has a la spread of 0.14 dex. Reddening is irrele- 
vant. 

• The index log[A^//6584]/[0//3727] is -0.99 ±0.07 
for Ay =0 and -1.19 ± 0.07 f or Ay = 0.7. The cal- 
ibration of Kcwlcv & Dopit^ (I2OO2D as modified by 
iKewlev fc EllisonI (2008) produces a double-valued 
abundance: 12 + \og{0/H) 8.21, 8.34 ± 0.06 for 
Ay = 0; and 12 + \og{0/H) ^ 8.21, 8.22 ± 0.06 for 
Ay = 0.7. 

• The Ne302 index, 
log([iVe///]A3869/[OJ/]A3727) = -0.96 ± 0.07, 
yields 12 + lo g(Q/g) = 8.1 9 ± 0.08 via the 
calibration of iShi et al] (|2007[ ). Reddening is 
irrelevant. This diagn ostic does not appear in 
IKewlev fc EllisonI ([20M ). so we cannot remove any 
calibration offsets. 

The R23 index, log[(A3727-|-A4959-|-A5007)/i?^], is un- 
fortunately double-valued as well. We measure logi?23 
as 0.94, 0'.96, 0.99 for A„=0.0, 0.7, 1.4, all with uncer- 
tainties of ±0.02 from the propagated flux uncertainty. 
(T his is similar to the value of 0.92 measured for cB58 
by iTephtz et all I2OOOI ) Such a high R23 requires a 
hi gh ionization parameter {\ogU > —2.87 from Fig. 5 
of IKewlev &: Dopital I2002D . The several R23 calibra- 
tions methods use various means to separate the "lower 
branch" an d "upper branch" , typi cally [N II] /[O II] or 
[N II] /Ha (jKewlev fc EllisonI I2OO8I ) . Unfortunately, for 
RCS0327 these ratios are on the border between upper 
and lower branch. Thus, we compute the abundance for 
each branch: 

• iZaritskv et al.l ()1994D published an R23 calibration 
for the upper branch only. We use it with caution 
since RCS0327 is not clearly in either the upper or 
lower branch. For Ay = 0, this calibration yields 
12 + \og{0/H) ^ 8.15 ± 0.02. For Ay = 6.7, it 
yields 8.14 ±0.02. 

• The iPilvugin fc ThuanI (|2005D method yields 12 + 
\og{0/H) = 8.26 ± 0.05 and 8.07 ± 0.08 for the 
upper and lower branches for Ay = 0, and 8.19 and 
8.22 (same uncertainties) for Ay = 0.7. 

• The iKobulnickv fc KewlevI ()2004[ ) method yields 
12±log(0/i?) = 8.19 and 8.16±0.02 for the upper 
and lower branches and Ay = 0, and 8. 17 ±0.02 for 
both branches at Ay = 0.7. 

These inferred abundances are plotted in Figure [SJ 

3.12. Abundances of other elements 

Following Ibotov et al.l (|2006[ ). we constrain the abun- 
dance of N+, Ne^+, Ar2+, and At^+ relative to H+. Un- 
like O, the full suite of ionization states is not observed 
for these elements, so we must apply ionization correc- 
tion factors to infer abundances. 

• N: We derive 12 + \og{N/H) > 6.53 (> 6.42) for 
Ay = 0.0 (0.7). Thus, N is depleted relative to 
the solar value by no more than 1.41 dex. The 
uncertainty from extinction is larger than from flux 
uncertainties. 
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• Ne: We derive 12+log(iVe/iJ) > 7.33±0.09, which 
is 0.60 dex below solar. 

• Ar: We derive 12 + log(Ar/iJ) > 5.54±0.18, which 
is 0.86 dex below solar. 

Though we have measured lower limits on all abun- 
dances, we can tightly constrain the key abundance ra- 
tios, which depend only weakly on T,,: 

• N/O: For the maximum allowed Te, log(A^/0) = 
-1.70 ± 0.02. If the A4363 flux is one-third the 
upper limit, then Tg is 70% of the maximum, and 
log(A^/0) = -1.89 ± 0.04. These ratios are 0.84 ± 
0.02 to 1.03 ± 0.04 dex below the solar abundance 
ratio. 

• Ne/O: For the maximum allowed Tg, log(iVe/0) = 
-0.89±0.09 for = 0, which is 0.13 dex below the 
solar ratio. For = 0.7 the ratio is log(A'^e/0) = 
-0.72±0.09, which is 0.04 dex above solar. At 70% 
of the maximum permitted Te, the constraints are 
from 0.03 dex below to 0.14 dex above solar. Thus, 
no matter the actual Tg, the Ne/O ratio must be 
solar-like. 

• Ar/O: For the maximum allowed Tg, \og{Ar/0) = 
—0.89 ± 0.18, which is 0.13 dex below the solar 
value. For 70% of the maximum permitted Te, 
\og{Ar/0) = -0.62 ±0.18, which is 0.14 dex above 
solar. 

Thus, RCS0327 has an abundance pattern in which 
the O abundance is at least 29% of solar, the Ne/O and 
Ar/O ratios are solar-like, and the N/O ratio is less than 
15% of the solar ratio. 

3.13. Joint constraints on physical conditions from 
photoionization models 

As a cross-check on these diagnostics, we run spec- 
tral synthesis and photoionization models. The UV 
spectra come from Starbu rst 99 (v5.1, web ver sion; 
iLeitherer et al.i il999 and Vazquez fc Leither'eill2005D . as- 
suming continuous star formation, an upper mass cutoff 
of 100 Mq, the default IMF parameters, the Padova AGB 
models, and with stellar metallicity set at 40% of solar. 
We feed these UV spectra to the photoionization code 
Cloudy version c08.00 (Ferland et al. 1998), running a 
grid of models of given electron density and starburst age, 
and in each model optimizing the metallicity and ion- 
ization parameter, as constrained by the measured line 
intensities relative to EE/?, assuming a foreground screen 
extinction Ay — 0.7, the measured uncertainties for 
lines contiguous with H/3, and 20% relative fluxing er- 
rors for noncontiguous lines. The abundance pattern is 
"H 11" region. Our grid of electron density covers ±1ct 
of the value measured in H'6M 

Cloudy computes an average electron temperature of 
Te = 1.19 to 1.24 X lO "^ K, which is consistent with the 
value inferred in ij3.10l Cloudy converges on an oxygen 
abundance in the range 12 + \og{0/H) = 8.07 to 8.13, 
and an ionization parameter in the range logf/ = —2.57 
to —2.42. Figure S] illustrates these constraints. Thus, 
Cloudy converges on an O abundance that is 0.1 dex 
lower than via the Te method, and a log U that is higher 



by 0.2 -0.3 dex than derived from the iKewlev fc Dopital 
(|200^ method. 

The source of this discrepancy may lie with differ- 
ent assumptions of abundance pattern, filling factor, 
atomic data, or ionizing stellar spectra compared to 
IKewlev fc Dopital (j2002 ). Resolving this moderate dis- 
crepancy is outside the scope of this paper at this time, 
but it provides a caution that one must be careful, when 
comparing derived physical properties of galaxies, to use 
the same methodology. In this paper, we will adopt the 
log U derived via the Kewl ev fc Dopital (|2002f ) method, 
so that we may compare to other lensed galaxies in the 
literature. 

Dust grains were not included in the Cloudy models. 
As a test, we added dust grains with the same metallicity 
and abundance pattern as the gas, which changed the 
inferred metallicity and ionization parameter by only 2% 
and 1%. 

4. DISCUSSION 

Since this is the first investigation of the rest-frame op- 
tical lines in RCS0327, only the brightest "knots" were 
targeted, representing only 32% ± 4% (random) ±1% 
(systematic) of the total light of the arc f ^33t . RCS0327 
was selected as the brightest and one of the largest lensed 
sources in a large survey; despite its exceptional appear- 
ance, SED and lensing modeling show it to be a fairly 
typical object (jWuvts et al.ll2010() . As such, it is likely 
that the high surface brightness portions of the source 
galaxy coincide with high magnification, and that this 
coincidence makes RCS0327 such a spectacular example 
of strong lensing. This supposition would explain the dif- 
fering star formation rates derived from 24/^m imaging of 
the main arc and counter-image, as noted in §3.7. More- 
over, because the lensed source s traddles the lensing tan- 
gential caustic (c.f. Figure 6 in lWuvts et al.|[2C)10r) . the 
lens model suggests that only some regions of the galaxy 
are represented in the giant arc, namely the core, and one 
of two apparent spiral arms, and may be highly magni- 
fied. Thus the giant arc is likely not a fair representation 
of the source galaxy as a whole. 

Since superb spectra can be quickly obtained for these 
bright knots, it is their physical conditions that we probe 
in this paper, with a caveat that these may not be "ordi- 
nary" regions of this galaxy. Once the arc has been fully 
mapped in these emission lines, and a high-resolution 
lensing map has been derived from HST imaging, we 
should understand how much the physical conditions 
vary across the source plane, and locate the brightest 
knots in the source plane. It will then be possible to fully 
contextualize our results in terms of the range of physical 
conditions across the spatial extent of RCS0327. 

For now, we assume that the knots for which we 
have spectra are representative. Since these are highest- 
quality rest-frame optical spectra ever obtained for a z=2 
galaxy, and yield precise measurements of the physical 
conditions, we now consider these measurements in the 
context of the literature, and the conditions under which 
stars form in the distant universe. 

4.1. Extinction 

The following E(B-V) extinctions have been mea- 
sured from B aimer decrements in lensed blue galaxies: 
0.27 for cB58 (Teplitz et al...200a) ; 0.28 ± 0.04 for the 
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Clone (iHainline et all [20091): .45 ± 0.04 for the Cos- 
mic Horseshoe (IHainline et all 2009); 0.59 ± 0.08 for 



J0900 ("Bi an et al 
If inkclstci n et al 



112010'); 0.67 ±0.21 for the 8 o'clock arc 
10009). Our measurement of E(B-V) 
= 0.23 ± 0.23 for RCS0327 is consistent with the lower 
part of this range, but the uncertainty is large, stem- 
ming mostly from the flux uncertainty in H7. A spec- 
trograph that obtains Ha and H/3 simultaneously, for ex- 
ample FIRE on Magellan or LUCIFER on LET, would 
provide a very precise measurement of the B aimer decre- 
ment in this galaxy, allowing a detailed comparison of 
the relative extinctions suffered by the gas and the stars. 

4.2. The reliability of O abundance diagnostics at z=2 

For RCS0327, the non-detection of [O HI] 4363 A sets 
a strict constraint on the oxygen abundance via the "di- 
rect" Te method: it must be no more than 0.55 dex (0.48 
dex) below the solar value for Ay = 0.7 (A„ = 0), or 28% 
and 33% as percentages of solar. Neon provides a cross- 
check on this abundance measurement. Since both Ne 
and O are alpha-group elements, their ratio should be 
solar-like, and indeed our measured ratio is solar with 
an uncertainty of ±0.14 dex. We believe this to be the 
first time the Ne/0 ratio has been measured at z ~ 2. 

We now accept the Tg-method's oxygen abundance 
measurement, and assess the reliability of the "bright 
line" diagnostics. 

These diagnostics are empirically calibrated at z=0 
against the method, or calibrated with photoioniza- 
tion models. As such, they implicitly assume densities 
and ionization parameters typical of nearby galaxies, 
which may not be appropriate for high-redshift galax- 
ies, for reasons of luminosity bias and evolution in the 
galaxy population. In addition, these diagnostics have 
significant relative offsets (up to 0.7 dex) at z=0 whose 
origins are not understood and which must be empirically 
calibrated away using SDSS galaxies ([Kewlev fc EUisoiJ 
120081) : it is not clear that for rapidly star-forming z=2 
galaxies these offsets will be the same. 

Thus, it is important to test these bright-line diag- 
nostics in situ at z=2. Figure [S] compares the O abun- 
dances inferred in i j3.11l for the bright-line diagnostics 
in RCS0327. For ease of comparison, we plot Ay=0 and 
Ay=0.7 separately. We divide the calibrations into five 
categories: [N II]/[0 II], Ne302, N2, 03N2, and R23 
methods. 

The [N II] /[O II] diagnostic (IKewlev fc DopitallMl 
is consistent with the Tg method, in part because its 
double- valued results and the large uncertainty covers 0.3 
dex of parameter space. Because this index has a large 
extinction correction, we cannot adequately test it until 
the extinction of RCS0327 is more precisely measured. 

The Ne302 diagnostic is consistent with the 
method, especially for Ay > 0. 

The N2 calibratio ns of iPettini fc Pagell ()2004D and 
iDenicolo et all (|200l are consistent with the lower limit 
from the Tg method, especially for Ay > 0. This test 
is especially importance, since this diagnostic is the one 
used to measure the evolution in the mass-metallicity 
relation from z=0 to z=2 (Erb et al. 2006) and z=3 
(|Mannucci et al.ll2009l: iMaiolino et al. 20081 

The 03N2 diagnostic ( jPettini fc Pagel 120041 ) yields an 
abundance that is 2.5cr below the lower limit for 
Ay = 0, and is thus inconsistent with the Tg limit. How- 



ever, moderate extinction {Ay = 0.7) makes the results 
consistent. The 03N2 result is lower than the abun- 
dances from the N2 and [N II] /[O II] indices. 

The R23 index has a large extinction correction. As 
such, it is difficult to assess the performance of this di- 
agnostic given the current uncertainty in the extinction 
of RCS0327. For Ay = only the upper branch of 
iPilvugin fc Thuaiil ([2005') i s consisten t with the re- 
sult; for Av = 0.7 the Pilvu gin fc Thu an (200|) method 
is also consistent. The iZaritskv et al.l ((l994 ) method pro- 
duces an inconsistent result, but this is unsurprising since 
it is only valid for high metallicity. 

This is the first test in situ at z=2 of the bright line 
abundance diagnostics for a star-forming galaxy of typi- 
cal metallicity. To summarize, the R23 and [N II] /[O II] 
methods are not inconsistent with the Te method so long 
as the extinction exceeds zero; the definitive test must 
await a more precise extinction measurement, since both 
diagnostics have large reddening corrections. The Ne302 
diagnostic is consistent with the Te method, though the 
errorbars are larger than for N2 and 03N2, since [Ne 
HI] is not a terribly bright line. The N2 and and 03N2 
methods are most constraining given the current extinc- 
tion measurement; of these, the N2 results are entirely 
consistent with the Te method, and the 03N2 diagnostic 
is consistent only given significant but plausible extinc- 
tion. 

The only other published Te detection at z=2 of which 
we are aware is in Lens22.3, a low-mass, low -metallicity 
galaxy behind AbeU 1689 (lYuan fc KewlevI [2009.) . For 
Lens22.3, the Te method yields 12±log(0/iJ) = 7.5±0.1 
for zero exti nction and 7.3 ± 0.1 f or high extinction. 
We apply the IPettini fc Pagell (|200l N 2 and Q3N2 cal- 
ibrations to the published line fluxes of I Yuan fc KewlevI 
(j2063) (N II is not detected). From N2 < -2.30 comes 
an abundance of 12 + \og{0/H) < 7.47 (7.59) using 
the third-order (linear) N2 calibration. The 03N2 in- 
dex, logC>3j y2 > 3.11, yield s 12 + log(C>/ij') < 7.73 on 
the original IPettini fc Pagell (|2004D scale, which is out- 
side the range that can be converted to the N2 frame 
(jKewlev fc Ellis on 2008\ These N2 and 03N2 measure- 
ments are consistent with the Te result. However, all 
these diagnostics are highly incon sistent with the R23 re - 
s ult of 12 ± log(C>/g) = 8.0-8.3 (lYuan fc Kewlevll2009l) . 

lYuan fc KewlevI (|2009f ) noted the discrepancy between 
the R23 and Te results, and attributed the fault to the Te 
method. Given the consistency between the N2, 03N2, 
and Te methods, we suggest that it may instead be R23 
that is failing in Lens22.3. 

4.3. The reliability of abundance diagnostics in other 
z=2 lensed galaxies 

Only in RCS0327 and Lens22.3 are the spectra suffi- 
cient to measure abundances via the Te method. How- 
ever, the bright line diagnostics have been used for a 
number of galaxies. A re-analysis of the mass-metallicity 
relation at z > 2 is outside the scope of this paper. How- 
ever, it is appropriate at this juncture to re-examine the 
N2 and 03N2 diagnostics in light of the results above. 
In Figure ini we plot the mass-metallicity relation at z=2, 
using literature results for lensed and unlensed galax- 
ies, and our results for RCS0327 . For fair comp arison, 
we convert stellar masses to the IChabrieil (l2003l) IMF, 
take line fiuxes from the literature and apply the N2 and 
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03N2 calibrations of iPettini fc Pagell ()2004l), brought to 
the N 2 system via the conversiorPI of IKewlev fc Elhsonl 
(|20n8h . 

We plot four le nsed galaxies in Figure IH cB58 
(|Teplitz et al |[2000l). J0900+ 2234 (B ian et al.|[20lo[ ). the 
8 o'clock arc (|Finkelstein et al... 2009 ). and RCS0327 (this 
work). These four galaxies span a range of 400 in stellar 
mass, ext ending below th e mass range probed by stacked 
samples ()Erb et alJ 120061 ). In aU four cases, the 03N2- 
derived abundances are systematically lower than those 
derived from N2, strikingly so in the case of the 8 o'clock 
arc. (Recall that the relative offset between these in- 
dices, as measured in SDSS at z=0, has already been 
removed.) These offsets were noted for each galaxy in 
their respective papers, but were dismissed as less than 
the calibration dispersion of the diagnostics. The Clone 
and Horseshoe lack stellar mass measurements to plot 
them on Figure IH but comparison of their N2 and 03N2- 
derived abundances shows the same trend: the 03N2 
index is low by 0.16 and 0.07 dex. 

It is clear from a consideration of all six lensed galaxies 
that the 03N2 and N2 indices are systematically offset 
at z=2. While this offset is small for the galaxies of 
intermediate metallicity, it is large for the 8 o'clock arc 
at the high metallicity end. 

This behavior is not surprising when we consider that 
03N2 has 5007/H/3 in the numerator, and that what 
makes z=2 galaxies stand out in the BPT diagram 
[Baldwin et"al1 ()1981[ ) is their high 5007/H/3 compared 
to SDSS galaxies. Thus, we suggest that high ionization 
conditions in z=2 galaxies cause the 03N2 diagnostic to 
function poorly at these redshifts. 

4.4. Lensed galaxies and the z^2 mass-metallicity 
relation 

RCS0327 has a stellar mass of log(MH.) = lO.OiO.l Mq 
(|Wuvts et al.ll2010[ ): both the main arc and the counter- 
arc give consistent results. This is 0.94 ± 0.14 dex lower 
than the Schec hter function parameter M*^^^ at z=2.0 
(|Marchesini et al. 2009 ), and lower than most of the Ly- 
man break galaxies in lErb et aTl ()2006D . The low mea- 
sured velocity dispersion qualitatively supports a low 
mass. In the mass-metallicity plane, RCS0327 lies sig- 
nificantly, 0.17 ± 0.04 dex, below the z=2 relation of 
lErb et aTl ()2006[ ) . This offset is too large to be measure- 
ment error, and so wc attribute it instead to a real abun- 
dance difference: RCS0327 is metal-poor by about 50% 
for its stellar mass and redshift, compared to lErb et al.l 
((20061) . 

J0 900-I-2234 shows this effect to an even greater de- 
gree (|Bian et al.ll2010 l) , since it has almost the same N2- 
derived oxygen abundance and a slightly higher stellar 
mass. 

On the high-mass end, the N2-derived oxygen abun- 
dance of the 8 o'c lock arc is quite consistent with the 
high-mass side of lErb et all (|2006D . perhaps 0.05-0.08 
dex low. By contrast, at the very low-mass end, cB8 
lies far off the relation, having a very low stellar mass 
but a high oxygen abundance, more typical of the mass- 
metallicity relation at z=0 than at z=2-3. Thus, even 
from a modest sample of four lensed galaxies, we are 



already beginning to probe the intrinsic spread in the 
mass-metallicity relation at z=2. 

4.5. Abundance pattern 



In ij3.12l we constrained the N/0 ratio as log(iV/0) < 
— 1.70. For plausible values of Tg the ratio can be lower 
by 0.2 dex. These abundances are startlingly close to 
the values measured for cB58: log(iV/0) = -1.76 ± 0.2 
dex IPc ttini ct al. 2002) and 12 + \og{0/H) = 8.26 
(jTepliTz et al.l 12000) F"! Even though RCS0327's stellar 
mass is 18 times larger than cB58's, and its redshift 
lower, star formation in these two galaxies has produced 
very similar rat ios of N, O, and H. UV spectroscopy for 
RCS0327, as in IPettini etall (|200l for cB58, should al- 
low element-by-element comparison of abundance ratios. 

For now, we concentrate on the N/0 ratio. RCS0327 
lies near the intersection of the primary and secondary N 
production lines. In other words, it lies right on the trend 
for secondary N production, and is somewhat below the 
primary plateau of \og{N/0) ^ —1.5. As such, its N/0 
ratio is at the low end of what has been measured for H 
II regions in spiral galaxies of comparable O ab undance, 
and is typical of dwarf galaxies (Ivan Zee et al.ll99a') . An- 
alytic models such as iHenrv et al.l ()2000h combine both 
N production mechanisms to produce smooth curves of 
N/0 versus O. RCS0327 lies on these curves near the 
"knee", where the N/0 ratio rapidly transitions from 
being independent of the oxygen abundance ( "primary" 
production), to being highly dependent on the O abun- 
dance ("secondary" production). C omparing RCS0327 
and cB58 to the numerical models of IHenrv et al.l ()2000f ) 
(their Figure 3b) suggests that both galaxies have rela- 
tively high star formation efficiencies. With a sample of 
two galaxies, it is premature to draw conclusions about 
how galaxies build up their nitrogen, but the strict N/0 
and O measurements we have made bode well for the 
future, if such work can be repeated for a larger sample. 

4.6. Ionization parameter 

Ionization parameters have been reported for a few 
other lensed galaxies, derived from the 5007/3727 di- 
agnostic. Unfortunately, this diagnostic has a strong 
dependence on O abundance, so this must be input to 
the diagnostic . We use the 50 07/3727 line fluxes re- 
ported by Hai nline et al.l ()2009f ) for the Horseshoe and 
the Clone, assume 40% solar metallicity (Asplund sys- 
te m), as inferred for both from the N2 index. Converting 
to lAnders &: Grevessd (Il989f) abundan ces and using the 
calibration of ' Kewlev fc Dopital (l200l . this corresponds 
to an ionization parameter of log U = —2.8 for the Horse- 
shoe and —2.7 for the Clone, without correcting for ex- 
tinction. Th is is consistent within 0. 1 dex with the values 
reported by iHainline et al.l ()2009f ). iTeplitz et all ()2000f ) 
did not calculate an ionization parameter for cB58, so 
we take thei r line fluxes and the O abundance from 
IPettini et all (|2002[ ). and find \ogU = -2.85 for cB58 
without correcting for extinction. Extinction dominates 
the uncertainty on this measurement, and acts to lower 
the ionization parameter. 

These ionization parameter measurements are entirely 
consistent with the value measured for RCS0327 of -2.73 



^'^ This conversion is modest (< 0.05 change in 12 + \og{0 / H)) 
at these metaUicities. 



" ITeplitz et aD |(2000D infer an N/O ratio t hat is higher by 0.5 
dex, a discrepancy discussed by IPettini et al.l 1)20021 ) . 
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to -2.85 (for = 0.7; 0.1 dex higher for A„ = 0). Thus, 
we conclude that the four lensed galaxies examined to 
date with this diagnostic all have very similar ionization 
parameters of about log U ^ —2.7. 

We now compare to local galaxies. The sample of 
iKewley et al.l (|2001[ ) would be ideal for comparison, since 
these are IR-luminous galaxies at Z'-~^0, but unfortunately 
their spectroscopy did not extend blueward to 3727 A. 
Instead, we turn to two papers which measured ioniza- 
tion p arameters for 65 SINGS galaxies (jMoustakas et al.l 
[Ml . 412 star-for ming galaxies, an d 120,000 galax- 
ies fr om the Sloan Digital Sky Survey (jMoustakas et al.l 
[20061) . using the same methodology of calculating log U 
as we do. The SINGS galaxies have a median and mode 
\ogU of —3. 0. The SPSS and 412 star-forming galaxies 
(figure 12 of iMoustakas et al1l2006[ ) have a narrow range 
of —3.2 < \ogU < —2.9. The prototypical starburst 
galaxy M82 has log U = -2.9 ± 0.07. 

Thus, the four lensed galaxies with measured ion- 
ization parameters all have the same measured value, 
within the uncertainties. This value, —2.7 to —2.8, is 
roughly twice as high as the median for samples of lower- 
luminosity nearby galaxies, and is 30-60% higher than 
for M82. Thus, the ionization parameters of these four 
Z'^2 star-forming galaxies are somewhat higher than lo- 
cal galaxies of much lower luminosity. 

4.7. Electron density 

In i j3.9l we derived an electron density of rig = 
252t^° cm-3 at 1.2 X 10* K from the [O II] 3726, 3729 
doublet ratio. This is currently the most precise such 
measurement at such high redshift. We now compare to 
density constraints from the literature for other lensed 
galaxies. 

.Hainline et al.l (poM ) measured the [S II] A6717, A6732 
flux ratio in two lensed galaxies. For the Clone, the 
result is 0.9 ±0.1 by extracting the fluxes in each of 
two apertures and summing, versus 0.75 ± 0.25 by ex- 
tracting all at once. For the Horseshoe, the result from 
the summed aperture is 1.0 ± 0.35. Using IRAF's tem- 
den at Te = lO'' K, these ratios indicate densities of 



900 



+500 
-300 



cm-'^ and 17001^^^^'' cm"^ for the Clone, and 
SOO^soo" cm-3 for ^Yie Horseshoe. iBian et al.l (pOlOl ) 
measured a [S II] flux ratio of 0.86 ± 0.2 for the sum 
of both apertures in J0900, which yields a density of 
llOOt^™" cm-3. In the rest-UV, IQuider et al.l (l2009l ) 
measure a C HI] flux ratio of /(1906)//(1908) = 1.1±0.2 
for the Horseshoe, corresponding to a density range of 
5000-25000 cm~^, which is inconsistent w ith the density 
range measured bv lHainline et al.l ()2009[ ). 

Thus, these literature measurements of electron den- 
sity are in the range 600-5000 cm~'^, albeit with large un- 
certainties. The low precision of these measurements in- 
dicates a clear need for deeper spectra to better measure 
electron density. Nevertheless, the current measurements 
for the Clone, Horseshoe, and J0900 favor high electron 
densities, much higher than the precise value we mea- 
sure for RCS0327. Thus, at present it is not clear what 
is a typical electron density for a star forming galaxy at 
these epochs. Additional high-quality measurements are 
urgently needed. 

4.8. Location in the BPT diagram 



The "BPT" diagnostic diagram of iBaldwin et al.l 
(|1981f) is commonly used to characterize the ionization 
conditions in galaxies. For RCS0327, we measure line 
ratios of log ([N II] 6583 / Ha) = -1.18 ± 0.07, and log 
(5007/H/3) = 0.69±0.02 for Av=0 and 0.13 less than that 
for = 0.7. In Figure [7] we plot RCS0327 on the BPT 
diagram; its high OIII/H/3 and extremely low NII/Ha 
place it in the upper left quadrant, close to the maximal 
starburst line. 

This is an exceptional position compared to the z=0 
SDSS, which has only 5 galaxies in that region of the 
BPT diagram. Ho wever, z=0 IR-lumin ous galaxies do 
occupy that space: iKewlev et al.l ()200lD have 9 galaxies 
with NII/Ha < —1. Since extreme star formation is 
rare in the local universe, these luminous star-forming 
galaxies may be a better basis for comparison than SDSS. 

It has previously been noted that z>l galaxies tend to 
be offs et toward higher 5007/H/3 ratios in t he BPT di- 
agram (IShaplev et al.ll2005l: lErb et al.l l2006l: [ Kriek et al.l 
120071) . iBrinchmann et al.l ()2008[) proposed that this is 
caused by an elevated ionization parameter at higher red- 
shift, and enumerated the following possible underlying 
causes: a top-heavy initial mass function; higher elec- 
tron densities; a higher volume filling factor; or a higher 
escape fraction of UV photons. 

Figure [7] shows that of the lensed galaxies, J0900 and 
the Clone show an offset similar to that of the lErb et aTl 
([2006) stacked galaxies, while t he 8 o'clock arc is consid- 
erably higher, as discussed bv Finkelstein et al.l (12009*) . 
By contrast, R CS0327 and cB58 are not as offset - they 
lie between the lErb et al.l ()2006[) points and the z=0 IR- 
luminous galaxies 

RCS0327 has a electron density (n^ = 235^26 cm-^) 
that is lower than the best-fltting densities for the Clone, 
Horseshoe, and J0900, though those other measurements 
have very large uncertainties. Its measured ionization pa- 
rameter (2.9 ± 0.17 for Ay = 0.7) is entirely consistent 
with measurements of the same diagnostic in the Horse- 
shoe, Clone, and cB58. Thu s, in two way s RCS 0327 con- 
tradicts the picture of Brinchmann et al.l ()2008D for BPT 
behavior at z=2: First, though its ionization parameter 
is the same as the other z=2 galaxies, it is not as offset 
in the BPT diagram; second, the electron density is not 
high, as has been suggested to explain offsets in the BPT 
diagram. 

Of the IR-luminous sample of IKewlev et al.l (|2001[ ) 
with NII/Ha < —1, the median electron density is 
~ 100 cm"'^. This is another case in which offsets in 
the BPT diagram do not appear to be caused by high 
density. 

At present, too few z=2 galaxies have the high-quality 
spectroscopy necessary to fully map their behavior in the 
BPT diagram, and link offsets back to evolution in phys- 
ical conditions. This will presumably change as new 
lensed galaxies are pursued, and as new multi-object 
near-IR spectrographs push down the luminosity func- 
tion of the non-lensed galaxy population. That said, 
RCS0327 demonstrates that offsets in the BPT diagram 
are not driven by higher ionization parameters and den- 
sities, as has been suggested. It also demonstrates that 

The Ho rseshoe is not plotte d because H/9 was contaminated 
by a skyhne HHainline et al 1120091 ). CB58 was not plotted because 
[S II] was not observed. 
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star-forming z=2 galaxies may have quite similar ioniza- 
tion parameters and densities to z=0 galaxies. 

Last in our discussion of the BPT diagram, we note 
that RCS0327's extreme location places it as far as pos- 
sible from the z=0 AGN locus. Thus, it is unlikely that 
the nebular emission of this galaxy is dominated by an 
AGN, tho ugh we have not ru led out a low-luminosity 
AGN (c.f. [G reene fc Holl2007t ). Upcoming observations 
with the Chandra X-ray Observatory should test this. 

5. CONCLUSIONS 

The spectra published here total 1.3 hr of integration 
time, divided into 20, 20, and 40 minute integrations 
over three spectral regions. These spectra demonstrate 
the power of gravitational lensing to explore the physical 
conditions of star formation at the epoch when most of 
the Universe's stars formed. 

For the second time in any galaxy at z ~ 2, and for 
the first in an average-metallicity galaxy, we tightly con- 
strain the O abundance using the "direct" Tf, method, 
thanks to the very constraining non-detection of [O 
III] A4363. We use this constraint to test the bright-line 
diagnostics of oxygen abundance, which are the easiest 
to measure, but are empirically calibrated at z=0 and 
thus incorporate assumptions about density and ioniza- 
tion parameter that may well be wrong at z~2. 

We find that the O abundance inferred from the N2 
index (the ratio of [N II] / Ha) agrees closely with the 
method and has a small uncertainty. Ne302 also per- 
forms well, albeit with a larger uncertainty since [Ne 
III] is not a bright line. N203 and R23 depend so 
strongly on the extinction that we cannot definitively as- 
sess their performance in RCS0327, though they appear 
to work for the current best estimate of the extinction, 
Al, ~ 0.7. T his is especially in teresting because R23 fails 
in Lens22.3 (|Yuan fc Ke wlcv 20091. the only published 
example of an [O III] 4363 detection at z^2. 

The 03N2 diagnostic is on the border of disagreeing 
with Tf. in RCS0327, depending on extinction. Compar- 
ing to the N2 index, 03N2 predicts systematically lower 
abundances in five Z'~2 lensed galaxies, dramatically so 
at near-solar metallicity. We suggest that 03N2 does 
not work in z=2 galaxies, perhaps due to different phys- 
ical conditions compared to z=0 where this diagnostic is 
calibrated. 

After ah, 03N2 is effectively a location in the BPT 
diagram, and it has been shown, via small samples of 
lensed galaxies and stacked samples of unlensed galaxies, 
that z=2 galaxies are offset in the BPT diagram, with 
higher 5007/H/3, than the cloud of SDSS galaxies at z=0. 
Indeed, RCS0327 is offset in the BPT diagram compared 
to the SDSS and iKewlev et al.l ()200lD galaxies, as offset 
as cB58 is, but not as offset as other z=2 galaxies. 

This is particularly interesting given that the measured 
ionization parameter for RCS0327 is the same as has been 
measured in the Horseshoe, the Clone, and cB58: log U ~ 
-2.7 to -2.8, which is 30-60% higher than in M82. 

We also tightly constrain the electron density: rig = 
252t^g cm-3 at 1.2 X lO'^ K, which is not terribly high 
compared to local galaxies. Thus, we conclude that it is 
premature to conclude that z ~ 2 star-forming galaxies 
have extremely high densities and ionization parameters, 
since the best-measured example, RCS0327, does not. 

We measure the relative abundances of N, Ne, and Ar 



compared to O. The Ne/0 and Ar/0 ratios are solar with 
uncertainties of ±0.14 dex, which is reassuring since all 
these elements are alpha-process and should enrich in 
lockstep. We believe this to be the first time the Ne/0 
ratio has been measured at z ~ 2. The N/0 ratio is one 
dex below the solar value, indicating that secondary N 
production has not yet begun in earnest. The O abun- 
dance and N/0 ratio are startlingly similar to those of 
cB58; it is not clear whether this agreement is merely 
coincidental, or indicates characteristic values for star- 
forming galaxies at this epoch. 

6. FUTURE DIRECTIONS 

This is by no means the last word on RCS0327 or on 
diagnostic spectroscopy of lensed galaxies. The follow- 
ing observations should significantly increase what can be 
learned about the physical conditions of RCS0327. First, 
a direct measurement [O III] 4363, or an even stricter 
upper limit, should provide a more stringent test of the 
bright-line O abundance diagnostics. This comparison, 
and a host of other constraints, are limited by the cur- 
rent uncertainty in the measured extinction. This would 
be improved by a deeper integration of H/3/II7 or a si- 
multaneous measurement of IIa/H/3. 

Thus far, we have considered only the spatially- 
integrated spectrum across the brightest portion of the 
arc. It will be fascinating to spatially map the physical 
conditions across this portion, and the fainter sections 
of the arc, to see how widely these physical parameters 
vary across the galaxy. Of course, such work requires a 
better lensing model, which will be enabled by pending 
HST observations. 

Finally, we humbly remember that a single galaxy can 
be a maverick, and that only by repeating this work in 
a representative sample of lensed galaxies will the phys- 
ical conditions of star formation at this epoch be confi- 
dently characterized. Larger samples will also fill in the 
BPT and mass-metallicity relations, to characterize the 
scatter in these relations and explore the reasons for the 
scatter. 
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Figure 2. The NIRSPEC spectra. The fluxed spectra are plotted in black, with the la error spectrum in red. The X-axis shows observed 
wavelength in Angstroms; the Y-axis shows observed flux in units of erg s~^ cm~^. Blue labels mark detected emission lines, and lines for 
which we measure upper limits. 



Physical conditions in a lensed z=1.7 galaxy 13 




1.005 1.010 1.015 
Observed Wavelength (micron) 

Figure 3. Two-component Gaussian fit to the [O II] 3726, 3729 doublet. Red lines show the fit to each emission line; a green line shows 
the summed fit. Blue shows the la error spectrum. 
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Figure 4. Starburst99 + Clo udy models, for constant star formation. A grid of models was run, varying the electron density through 
itlcr of the value measured in i|3.9l and the starburst age. The input spectrum changes slightly due to stochastic effects. The filled circles 
show models with the density measured assuming Tg = 1.2 X 10^ K, and open circles show models using the density measured assuming 
Te = 10^ K. The measured spectral lines and the density measurement tightly constrain the ionization parameter and metallicity. 
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Figure 5. Comparison of oxygen abundance diagnostics for RCS0327. The top pan el assumes = 0, and the bottom panel assumes 
Ay = 0.7. The relative calibration offsets observed at z=0 IIKewlev fc EllisonI 120081) have already been removed; all indices are on the 
relative frame of N2 in Pettini fc Pagel (12004 ). 
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Figure 6. The Mass-Metallicity relation at z=0 and z=2. SDSS galaxies define the z=0 relation (grey cloud). The z=2 relation is shown 
by sta cked spect ra of unlensed z=2 LBGs from Erb ct al. (2006) (hollow squares), as well as four lensed galaxi es at z~2: J0900+2234 
from Bian et al.l l|20ia ): cB58 from iTeplitz et a l. (2000) and Siana et al. (2008); the 8 o'clock arc from Fi nkelstein et al.l p00<3); and 
RC S0327 fr om this work (filled points). For consistency, when authors fit stellar masses using a Salpeter IMF, we have converted to a 
IChabricr ( 2003) IMF by dividing M« by 1.8. Squares and circles show abundances determined from the N2/Ha and 03N2 calibrations 
of Pcttini & Pagel (2004), respectively, where the 03N2 abundances have been brought onto the N2/Ho relative system using the small 
conversion of Kcwlev Ellison 12008 ). An arrow shows the lower limit on abundance derived from the T^ method for RCS 0327. Errorbars 
on the data show the propagated flux uncertainty. The errorbars at upper left show the Icr calibration uncertainties of IPettini &: PageH 
((2001). 
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Figure 7. The BPT diagram. The z=0 relation is defined by the SDSS DR7 spectroscopic sample dAbazaiian et ani2009l; black dots) , 
and by the IR-bright galaxies of Kcwlcy ct al. (2001) (orange a sterisks). Sta cked sampl es of galaxies at z=2 are from lErb et al.l ( 12009!) 
(blue diamonds); lensed z=2 galaxies are from .Bian et al.. t2010il . ri'eplitz et al.l ((2000i) . and ll^'inkelstein et al.l ( 1200911 (green diamonds), and 
this work (red diamond). 



